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Abstract 

The T2K experiment observes indications of — > v e appearance in data accumulated with 
1.43 x 10 20 protons on target. Six events pass all selection criteria at the far detector. In a three- 
flavor neutrino oscillation scenario with \ Amj 3 \ = 2.4 x 1(T 3 eV 2 , sin 2 2 6> 23 = 1 and sin 2 29 13 = 0, 
the expected number of such events is 1.5±0.3(syst.). Under this hypothesis, the probability to 
observe six or more candidate events is 7xlCT 3 , equivalent to 2.5a significance. At 90% C.L., 
the data are consistent with 0.03(0.04)< sin 2 29 13 < 0.28(0.34) for Sqp = and normal (inverted) 
hierarchy. 

PACS numbers: 14.60.Pq,13.15.+g,25.30.Pt,95.55.Vj 
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We report results of a search for v e appearance in the T2K experiment pQ. In a three- 
neutrino mixing scenario, flavor oscillations are described by the PMNS matrix [21 E] , usually 
parametrized by the three angles #i 2 , #23, #13, and the CP- violating phase #cp- Previous 
experiments have observed neutrino oscillations driven by #12 and #23 in the solar (Am^) 
and atmospheric (Am^ ^Aml) sectors [IH9]. In the atmospheric sector, data are consistent 
with |Am2 3 | — 2.4 x 10~ 3 eV 2 , a normal Am^ > or inverted Am^ < mass hierarchy, 
and sin 2 26* 2 3 close to, or equal to unity. Searches for oscillations driven by 6^3 have been 
inconclusive and upper limits have been derived [TUHT3"] . with the most stringent being 
sin 2 2fl 13 <0.15 (90%C.L.), set by CHOOZ [2] and MINOS [IS]. 

T2K uses a conventional neutrino beam produced at J-PARC and directed 2.5° off-axis to 
Super-Kamiokande (SK) at a distance L = 295 km. This configuration produces a narrow- 
band beam [16], tuned at the first oscillation maximum E u = |Am2 3 | L/(2ir) ~ 0.6 GeV, 
reducing backgrounds from higher energy neutrino interactions. 

Details of the T2K experimental setup are described elsewhere [17] . Here we briefly review 
the components relevant for the v e search. The J-PARC Main Ring (MR) accelerator [18] 
provides 30 GeV protons with a cycle of 0.3 Hz. Eight bunches are single-turn extracted in 
5 /j,s and transported through an extraction line arc defined by superconducting combined- 
function magnets to the production target. The primary beamline is equipped with 21 
electrostatic beam position monitors (ESM), 19 segmented secondary emission monitors 
(SSEM), one optical transition radiation monitor (OTR) and five current transformers. The 
secondary beamline, filled with He at atmospheric pressure, is composed of the target, 
focusing horns and decay tunnel. The graphite target is 2.6 cm in diameter and 90 cm 
(1.9\ int ) long. Charged particles exiting the target are sign selected and focused into the 
96 m long decay tunnel by three magnetic horns pulsed at 250 kA. Neutrinos are primarily 
produced in the decays of charged pions and kaons. A beam dump is located at the end of 
the tunnel and is followed by muon monitors. 

The Near Detector complex [T7j located 280 m downstream from the target hosts two 
detectors. The on-axis Interactive Neutrino GRID (INGRID) accumulates neutrino inter- 
actions with high statistics to monitor the beam intensity, direction and profile. It consists 
of 14 identical 7-ton iron-absorber/scintillator-tracker sandwich modules arranged in 10 m 
by 10 m crossed horizontal and vertical arrays centered on the beam. The off-axis detector 
reconstructs exclusive final states to study neutrino interactions and beam properties corre- 
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sponding to those expected at the far detector. Embedded in the refurbished UAl magnet 
(0.2 T), it consists of three large volume time projection chambers (TPCs) [19 j interleaved 
with two fine-grained tracking detectors (FGDs, each 1 ton), a 7r°-optimized detector and a 
surrounding electromagnetic calorimeter. The magnet yoke is instrumented as a side muon 
range detector. 

The SK water Cherenkov far detector [20] has a fiducial volume (FV) of 22.5 kton within 
its cylindrical inner detector (ID). Enclosing the ID all around is the 2 m-wide outer detector 
(OD). The front-end readout electronics allow for a zero-deadtime software trigger. Spill 
timing information, synchronized by the Global Positioning System (GPS) with < 150 ns 
precision, is transferred online to SK and triggers the recording of photomultiplier hits within 
±500 fis of the expected arrival time of the neutrinos. 

The results presented in this Letter are based on the first two physics runs: Run 1 (Jan- 
Jun 2010) and Run 2 (Nov 2010-Mar 2011). During this time period, the MR proton beam 
power was continually increased and reached 145 kW with 9 x 10 13 protons per pulse. The 
targeting efficiency was monitored by the ESM, SSEM and OTR and found to be stable 
at over 99%. The muon monitors provided additional spill-by-spill steering information. A 
total of 2,474,419 spills were retained for analysis after beam and far detector quality cuts, 
yielding 1.43 x 10 20 protons on target (p.o.t.). 

We present the study of events in the far detector with only a single electron-like (e- 
like) ring. The analysis produces a sample enhanced in v e charged-current quasi-elastic 
interactions (CCQE) arising from — > v e oscillations. The main backgrounds are intrinsic 
v e contamination in the beam and neutral current (NC) interactions with a misidentified 
7T°. The selection criteria for this analysis were fixed from Monte Carlo (MC) studies before 
the data were collected, optimized for the initial running conditions. The observed number 
of events is compared to expectations based on neutrino flux and cross-section predictions 
for signal and all sources of backgrounds, which are corrected using an inclusive charged- 
current (CC) measurement in the off-axis near detector. 

We compute the neutrino beam fluxes (Fig. [T]) starting from models and tuning them to 
experimental data. Pion production in (p, 9) bins is based on the NA61 measurements [2T] , 
typically with 5-10% uncertainties. Pions produced outside the experimentally measured 
phase space, as well as kaons, are modeled using FLUKA [221 [23]. These pions are assigned 
systematic uncertainties on their production of 50%, while kaon production uncertainties, 
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FIG. 1. Predicted neutrino fluxes at the far detector, in absence of oscillations. The shaded boxes 
indicate the total systematic uncertainties for each energy bin. 

estimated from a comparison with data from Eichten et al. [24], range from 15% to 100% 
depending on the bin. GEANT3 [25], with GCALOR ]26j for hadronic interactions, handles 
particle propagation through the magnetic horns, target hall, decay volume and beam dump. 
Additional errors to the neutrino fluxes are included for the proton beam uncertainties, 
secondary beamline component alignment uncertainties, and the beam direction uncertainty. 

The neutrino beam profile and its absolute rate (1.5 events/10 14 p.o.t.) as measured by 
INGRID were stable and consistent with expectations. The beam profile center (Fig. [2]) 
indicates that beam steering was better than ±1 mrad. The correlated systematic error is 
±0.33(0.37) mrad for the horizontal(vertical) direction. The error on the SK position relative 
to the beamline elements was obtained from a dedicated GPS survey and is negligible. As 
shown in Fig. [lj the estimated uncertainties of the intrinsic and v e fluxes below 1 GeV 
are around 14%. Above 1 GeV, the intrinsic v e flux error is dominated by the uncertainty 
on the kaon production rate with resulting errors of 20-50%. 

The NEUT MC event generator [28], which has been tuned with recent neutrino inter- 
action data in an energy region compatible with T2K [29H3T] . is used to simulate neutrino 
interactions in the near and far detectors. The GENIE [32J generator provides a separate 
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FIG. 2. Beam centering stability in horizontal (x, South-North) and vertical (y, Down-Up) direc- 
tions as a function of time, as measured by INGRID. Errors shown are only statistical. The dashed 
lines correspond to a change of beam direction by ±1 mrad. 

TABLE I. Summary of systematic uncertainties for the relative rate of different charged-current 
(CC) and neutral-current (NC) reactions to the rate for CCQE. 



Process Systematic error 

CCQE energy-dependent (7% at 500 MeV) 

CC lvr 30% (E v < 2 GeV) - 20% (E v > 2 GeV) 

CC coherent vr 1 * 1 100% (upper limit from [27] ) 

CC other 30% {E v < 2 GeV) - 25% (E u > 2 GeV) 

NC lvr 30% {E v < 1 GeV) - 20% (E u > 1 GeV) 

NC coherent vr 30% 

NC other vr 30% 

FSI energy-dependent (10% at 500 MeV) 



cross-check of the assumed cross-sections and uncertainties, and yields consistent results. A 
list of reactions and their uncertainties relative to the CCQE total cross-section is shown in 
Table [[} An energy-dependent error on CCQE is assigned to account for the uncertainty in 
the low energy cross-section, especially for the different target materials between the near 
and far detectors. Uncertainties in intranuclear final state interactions (FSI), implemented 
with a microscopic cascade model [33], introduce an additional error in the rates (see e.g. 
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[35). 

An inclusive CC measurement in the off-axis near detector is used to constrain the 
expected event rate at the far detector. From a data sample collected in Run 1 and corre- 
sponding to 2.88 x 10 19 p.o.t. after detector quality cuts, neutrino interactions are selected 
in the FGDs with tracks entering the downstream TPC. The most energetic negative track 
in the TPC is chosen and we require its ionization loss to be compatible with a muon. To 
reduce background from interactions outside the FGDs, there must be no track in the up- 
stream TPC. The analysis selects 1529 data events (38% CC efficiency for 90% purity, 
estimated from MC). The momentum distribution of the selected muons (Fig. [3]) shows good 
agreement between data and MC. The measured data/MC ratio is 

R^ ata /RND C = 1-036 ± 0.028(stat.)+°;^(det.syst.) ± 0.038(phys.syst.), (1) 

where i?^ ata and Rne> C are the P-o.t. normalized rates of CC interactions in data and 
MC. The detector systematic errors mainly come from tracking and particle identification 
efficiencies, and physics uncertainties are related to the interaction modeling. Uncertainties 
that effectively cancel between near and far detectors were omitted. 

At the far detector, we extract a fully-contained fiducial volume (FCFV) sample by 
requiring no event activity in either the OD or in the 100 fis before the event trigger time, at 
least 30 MeV electron-equivalent energy deposited in the ID (defined as visible energy E vis ), 
and the reconstructed vertex in the fiducial region. The data have 88 such FCFV events 
that are within the timing range from —2 to 10 /is around the beam trigger time. The 
accidental contamination from non-beam related events is determined from the sidebands 
to be 0.003 events. A Kolmogorov-Smirnov (KS) test of the observed number of FCFV 
events as a function of accumulated p.o.t. is compatible with the normalized event rate 
being constant (p-value=0.32). The analysis relies on the well-established reconstruction 
techniques developed for other data samples [I]. Forty-one events are reconstructed with 
a single ring, and eight of those are e-like. Six of these events have E vis > 100 MeV 
and no delayed-electron signal. To suppress misidentified ir° mesons, the reconstruction of 
two rings is forced by comparison of the observed and expected light patterns calculated 
under the assumption of two showers [35J, and a cut on the two-ring invariant mass Mj n „ < 
105 MeV/c 2 is imposed. No events are rejected (Fig. |4|. Finally, the neutrino energy 
El ec is computed using the reconstructed momentum and direction of the ring, by assuming 
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FIG. 3. Measured muon momentum of z/„ CC candidates reconstructed in the FGD target. The 
data are shown using points with error bars (statistical only) and the MC predictions are in 
histograms shaded according to their type. 

quasi-elastic kinematics and neglecting Fermi motion. No events are rejected by requiring 
E^ ec < 1250 MeV, aimed at suppressing events from the intrinsic v e component arising 
primarily from kaon decays (Fig. [5]). The data and MC reductions after each selection 
criterion are shown in Table [TTJ The v e appearance signal efficiency is estimated from MC 
to be 66% while rejection for + CC, intrinsic v e CC, and NC are > 99%, 77%, and 
99%, respectively. Of the surviving background NC interactions constitute 46%, of which 
74% are due to 7r° mesons and 6% originate from single gamma production. 

Examination of the six data events shows properties consistent with v e CC interactions. 
The distribution of the cosine of the opening angle between the ring and the incoming beam 
direction is consistent with CCQE events. The event vertices in cylindrical coordinates 
(R,(f),z) show that these events are clustered at large R, near the edge of the FV in the 
upstream beam direction. A KS test on the R 2 distribution of our final events yields a 
p-value of 0.03. If this was related to contamination from penetrating particles produced 
in upstream neutrino interactions, then the ID region outside the FV should show evidence 
for such events, however this is not observed. In addition, an analysis of the neutrino 
interactions occurring in the OD volume is consistent with expectations. 
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FIG. 4. Distribution of invariant mass Mj ni) when each event is forced to be reconstructed into two 
rings. The data are shown using points with error bars (statistical only) and the MC predictions 
are in shaded histograms, corresponding to oscillated v e CC signal and various background sources 
for sin 2 2#i3 = 0.1. The last bin shows overflow entries. The vertical line shows the applied cut at 
105 MeV/c 2 . 

To compute the expected number of events at the far detector N^g, we use the near 
detector CC interaction rate measurement as normalization, and the ratio of expected 
events in the near and far detectors, where common systematic errors cancel. Using Eq. [T], 
this can be expressed as: 



ly SK — \ n ND f n ND J ' 1V SK 5 



(2) 



where N^f is the MC number of events expected in the far detector. Due to the correlation 
of systematic errors in the near and far detector samples, Eq.[2]reduces the uncertainty on the 
expected number of events. Event rates are computed incorporating three-flavor oscillation 
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FIG. 5. Same as Fig. [4] for the reconstructed neutrino energy spectrum of the events which pass 
all v e appearance signal selection criteria with the exception of the energy cut. The vertical line 
shows the applied cut at 1250 MeV. 



probabilities and matter effects with Am 2 2 = 7.6 x 10 5 eV 2 , Am^ = +2.4 x 10 3 eV 2 , 
sin 26 12 = 0.8704, sin 2 2^3 = 1.0, an average Earth density p=3.2 g/cm 3 and Sep = unless 
otherwise noted. The expectations are 0.03(0.03) + CC, 0.8(0.7) intrinsic v e CC, and 
0.1(4.1) z/ M — > v e oscillation events for sin 2 2#i3=0(0.1), and 0.6 NC events. As shown in Ta- 



ble III, the total systematic uncertainty on N e ^ depends on #13. Neutrino flux uncertainties 
contribute 14.9%(15.4%) to the far(near) event rates, but their ratio has an 8.5% error due to 
cancellations. The near detector CC selection efficiency uncertainty yields 1J*2% and the 
statistical uncertainty gives 2.7%. The errors from cross-section modeling are dominated 
by FSI uncertainties and by the knowledge of the a"(z/ e )/a"(t' At ) ratio, estimated to ±6%. 
The systematic uncertainties due to event selection in SK were studied with cosmic-ray 
muons, electrons from muon decays, and atmospheric neutrino events. Their contribution 
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TABLE II. Event reduction for the v e appearance search at the far detector. After each selection 
criterion is applied, the numbers of observed (Data) and MC expected events of v„ CC, intrinsic 
v e CC, NC, and the v e CC signal, are given. All MC CC samples include three-flavor oscillations 
for sin 2 26>i 3 =0.1 and 5 C p = 0. 

Data z^CC ^ e CC NC ^v e CC 



(0) interaction in FV 


n/a 


67.2 


3.1 


71.0 


6.2 


(1) fully-contained FV 


88 


52.4 


2.9 


18.3 


6.0 


(2) single ring 


41 


30.8 


1.8 


5.7 


5.2 


(3) e-like 


8 


1.0 


1.8 


3.7 


5.2 


(4) E vis > 100 MeV 


7 


0.7 


1.8 


3.2 


5.1 


(5) no delayed electron 


6 


0.1 


1.5 


2.8 


4.6 


(6) non-7r°-like 


6 


0.04 


1.1 


0.8 


4.2 


(7) E' v ec < 1250 MeV 


6 


0.03 


0.7 


0.6 


4.1 



to SNg X ^/Ng X ^ for e.g. sin 2 2# 13 = 0.1 is as follows: 1.4% from the fiducial volume definition, 
0.6% from the energy scale and 0.2% from the delayed electron signal tagging efficiency. The 
7t° rejection efficiency, studied with a NC 7r° topological control sample combining one data 
electron and one simulated gamma event, contributes 0.9%. The uncertainty on the accep- 
tance of one-ring e-like events was studied with an atmospheric neutrino sample, adding a 
contribution of 5% from ring counting and 4.9% from particle identification uncertainties. 
The performance of muon rejection by the ring particle identification algorithm was inves- 
tigated using cosmic-ray muons and atmospheric neutrino events, giving 0.3%. The effect 
from uncertainties in the Mi nv cut is 6.0%. Combining the above uncertainties, the total far 
detector systematic error contribution to SNt^ /N^ is 14.7%(9.4%) for sin 2 26*i3 = 0(0.1). 

Our oscillation result is based entirely on comparing the number of v e candidate events 
with predictions, varying sin 2 29i 3 for each <5cp value. Including systematic uncertainties, 
the expectation is 1.5±0.3(5.5±1.0) events for sin 2 2^13 = 0(0.1). At each oscillation pa- 
rameter point, a probability distribution for the expected number of events is constructed, 
incorporating systematic errors [37], which is used to make the confidence interval (Fig. [6]), 
following the unified ordering prescription of Feldman and Cousins [38J. 

In conclusion, the observation of six single ring e-like events exceeds the expectation of 
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TABLE III. Contributions from various sources and the total relative uncertainty for sin 2 2#i 3 =0 
and 0.1, and Sep = 0. 





sin 2 2#i3 = sin 2 2% = 0.1 


(1) neutrino flux 


± 8.5% 


± 8.5% 


(2) near detector 


+5.6o/ 
-5.2/0 


+5.60/ 
-5.2/0 


(3) near det. statistics 


± 2.7% 


± 2.7% 


(4) cross section 


± 14.0% 


± 10.5% 


(5) far detector 


± 14.7% 


± 9.4% 


Total 8N e *g/N e *g 


+22. 80/ 
-22.7 /0 


+17.60/ 
-17.5/0 



a three-flavor neutrino oscillation scenario with sin 2 2#i 3 = 0. Under this hypothesis, the 
probability to observe six or more candidate events is 7xlCT 3 . Thus, we conclude that 
our data indicate v e appearance from a neutrino beam. This result converted into a 
confidence interval yields 0.03(0.04) < sin 2 20 13 < 0.28(0.34) at 90% C.L. for sin 2 20 23 = 1.0, 
|Am| 3 | = 2.4 xlO -3 eV 2 , <5cp = and for normal (inverted) neutrino mass hierarchy. Under 
the same assumptions, the best fit points are 0.11(0.14), respectively. For non-maximal 
sin 2 2#23, the confidence intervals remain unchanged to first order by replacing sin 2 2#i 3 by 
2 sin 2 #23 sin 2 26> 13 . More data are required to firmly establish v e appearance and to better 
determine the angle 9i 3 . 
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sin 2 26 13 

FIG. 6. The 68% and 90% C.L. regions for sin 2 26^ for each value of 5qp, consistent with 
the observed number of events in the three-flavor oscillation case for normal (top) and inverted 
(bottom) mass hierarchy. The other oscillation parameters are fixed (see text). The best fit values 
are shown with solid lines. 
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